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We show that axinos produced in the early Universe in the 
decay of the lightest neutralinos are a natural candidate for 
cold dark matter. We argue that axinos may well provide the 
main component of the missing mass in the Universe because 
their relic density is often comparable with the critical density. 

PACS: 14.80.-j, 14.80.Ly, 95.35 



1. Introduction. Axinos are predicted to exist in mod- 
els involving low-energy supersymmetry (SUSY) and 
the Peccei-Quinn solution M to the strong CP prob- 
lem. They are supersymmetric partners of axions 
SUSY is widely considered as perhaps the most attractive 
framework in which the Fermi scale can be naturally con- 
nected with physics around the Planck scale. The Peccei- 
Quinn (PQ) mechanism, which invokes a global, chiral 
U(l) symmetry group spontaneously broken at some high 
energy scale f a ~ 10 11 GeV remains the most compelling 
way of solving the strong CP problem. 

Axinos are thus very strongly motivated. Despite this, 
they have received much less attention in the literature 
than other SUSY partners. Of particular importance to 
both experimental searches and cosmology is the light- 
est supersymmetric particle (LSP). Axinos, being mas- 
sive and electrically and color neutral are an interesting 
candidate for the LSP. One of the most important conse- 
quences of supersymmetry for cosmology in the presence 
of unbroken R-parity is the fact that the LSP is stable 
and may contribute substantially to the relic mass den- 
sity in the Universe. If the contribution is of order the 
critical density p cr u, such a particle is considered an at- 
tractive dark matter (DM) candidate. Current models 
of the formation of large structures as well as measured 
shape of their power spectrum strongly suggest that a 
dominant contribution to the dynamical component of 
the total mass-energy density is that from cold DM ||. 

In the minimal SUSY model (MSSM), the LSP is usu- 
ally assumed to be the lightest of the four neutralinos. 
The (lightest) neutralino % is a mixture x — Z\\B + 
Z12W3 + + Zi^II® of the respective fcrmionic part- 

ners (denoted by a tilde) of the electrically neutral gauge 
bosons B and W3, and Higgs bosons Hi, and H u . It is 
well-known that the neutralino's relic density p x is of- 
ten of order p cr it- At high temperatures in the early 
Universe a thermal population of neutralinos remains in 



equilibrium with the thermal bath. When their annihi- 
lation rate into ordinary matter becomes smaller than 
the expansion of the Universe, they decouple from the 
thermal bath, or "freeze-out" j5j. The "freeze-out" tem- 
perature is typically small compared to the neutralino's 
mass, Tf ~ m x /20 ||. Relic neutralinos are therefore 
always non-relativistic, or cold, DM candidate 

Many bounds on the neutralino mass and the param- 
eter space of SUSY models have been derived by re- 
quiring that the neutralino abundance does not "over- 
close" the Universe. This requires satisfying the condi- 
tion Vt x h 2 < 1 where £l x = p x /p cr it and h is related 
to the Hubble parameter Hq = 100 /ikm/sec/Mpc. This 
condition comes from considering the evolution of a ther- 
mal population of LSPs in the expanding Universe and in 
particular their annihilation cross section at decoupling. 
The annihilation has to be efficient enough to deplete the 
LSP number density to acceptable values. Cosmological 
properties of the neutralino as the LSP and DM are often 
taken into account in many studies of SUSY, including 
present and future collider and DM searches. 

In this Letter we will show that this standard paradigm 
changes dramatically if one assumes that it is the axino, 
rather than the lightest neutralino, which is the LSP. 
This assumption is well justified. Experimental searches 
at LEP have now pushed the neutralino mass limit con- 
siderably, above about 28 GeV in the MSSM Q. In 
more restrictive, and perhaps more motivated, models 
the bound can be much higher. For example, in the Con- 
strained MSSM (CMSSM) [§, also known as the minimal 
supergravity model, it is already around 42 GeV |jl|| . On 
the other hand, it is worth remarking that these bounds 
strongly depend on the (well-motivated) assumption that 
the masses of the gauginos (the fermionic partners of 
the gauge bosons) are equal at a grand-unified scale. 
In the absence of this condition one recovers a robust 



model independent bound m x > 3 GeV 10 from requir- 



ing il x h 2 < 1 - a neutralino version of the so-called Lee- 
Weinberg bound ||. 

In contrast to the neutralino, the mass of the axino, 
m~, remains not only virtually unconstrained experimen- 
tally but also theoretically one can easily imagine it in 
the few to tens of GeV range which we are interested 
in |[^]. This is illustrated by the following examples. 
In the supersymmetric version of the heavy quark axion 
model (the KSVZ model |]), the axino mass can arise 
at one-loop level with a SUSY breaking ^4-term insertion 
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at the intermediate heavy squark line. Then, we expect 
to- ~ {Jq/Stt 2 )A where /q is the Yukawa coupling of the 
heavy quark to a singlet field containing the axion. In 
a straightforward SUSY version of the DFSZ || model 
axino mass is typically rather small, to- ~ O(kcV) as 
has been pointed out in Ref. ]l3| ]. However, in addi- 
tion to the above inevitable contributions to the KSVZ 
and DFSZ axino masses, there can be other contributions 
from superpotentials involving singlet fields. For exam- 
ple, if the PQ symmetry is assumed to be broken by 
the renormalizable superpotential term (KSVZ or DFSZ 
models) W = fZ(SiS% — f 2 ) where / is a coupling, and 
Z, Si and S 2 are chiral fields with PQ charges of 0, +1 
and —1, respectively, then the axino mass can be at the 
soft SUSY breaking mass scale. The axino mass arises 
from the mass matrix of Si, S 2 , and Z, 



0, m~, ffa 
m~, 0, ff a 

ffa, ffa, 



(1) 



Certainly, the tree level axino mass is zero if < Z >= 0. 
However, with soft terms included, there appears a linear 
term in Z, V = \f\ 2 (\Si\ 2 + \S 2 \ 2 )\Z\ 2 + {AfSiS 2 Z+h.c); 
thus < Z > is of order A/f, and the axino mass can 
arise at the soft mass scale. A complete knowledge of 
the superpotential is necessary to pin down the axino 
mass |l^,|l4|]. Therefore, generically it is not unreasonable 
to consider the axino mass scale of order tens of GeV. 

One severe bound of to— < 2 keV has been derived by 
Rajagopal, Turner and Wilczek |l3|]. This bound arises 
from requiring that the "primordial" axinos, produced 
along with the axions in the very early Universe when 
the PQ symmetry becomes broken around the scale f a ~ 
10 11 GeV, do not contribute too much to the total relic 
abundance of the Universe, Q~h 2 < 1. As was noted in 
Rcf. jl3) , the bound to— < 2 keV (which would make the 
axino a warm dark matter candidate) can be evaded by 
assuming that, at temperatures below f a , the Universe 
underwent a period of inflation and that the temperature 
of subsequent reheating was sufficiently below f a . These 
assumptions are not too radical and have now become 
part of the standard cosmological lore || . The same way 
out is also usually used to solve the analogous problem 
with primordial gravitinos Jig] . 

Once the number density of the primordial axinos has 
been diluted by inflation, they can be again produced in 
the decays of heavier particles Since axino's cou- 

plings to matter are strongly suppressed by 1/ f a , all 
heavier SUSY partners first cascade-decay to the next- 
to-lightest SUSY partner (NLSP). A natural candidate 
for the NLSP is the lightest neutralino. As stated above, 
the neutralino "freezes-out" at Tf ~ to x /20. If it were 
the LSP, its co-moving number density after freeze-out 
would remain nearly constant. In the scenario considered 
here, the neutralino, after decoupling from the thermal 



equilibrium, will subsequently decay into the axino via, 
e.g., the process 



(2) 



This process was already considered early on in 
Ref. @ (see also fT|) in the limit of a photino NLSP 
and only for both the photino and axino masses assumed 
to be very low, to— < 1 GeV and to— < 300 eV, the for- 
mer bound now excluded by LEP searches. In that case, 
the photino lifetime was typically much larger than 1 
second thus normally causing destruction of primordial 
deuterium produced during nucleosynthesis by the ener- 
getic photon. Avoiding this led to a lower / a -dependent 
bound on the mass of the photino in the MeV range [^6| . 

In this Letter, we show that the conclusions and 
bounds of Refs. Jl(],[L| can be evaded if one considers 
both the axino and the neutralino in the GeV mass 
range. In this regime the neutralino decays into the ax- 
ino typically well before nucleosynthesis thus avoiding the 
problems considered in Refs. The resulting non- 

thermally produced axino will be a cold DM candidate. 

2. NLSP Freeze- Out. The effective coupling of the 
neutralino with the axino is very much weaker than 
that of its interactions with other matter. Therefore 
the neutralino's decoupling is not different from the case 
when it is the LSP. The freeze-out temperature Tf is 
determined by the annihilation cross section o~(xx ~ * 
ordinary matter) and is normally well-approximated by 
iteratively solving the equation for Xf — Tf /m x 
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where M P = 1.22 x 10 19 GeV, N F is the effective num- 
ber of relativistic degrees of freedom and (o~v re i) is the 
averaged product of annihilation cross section and the 
annihilating neutralinos' relative velocity. As already 
mentioned, the scaled freeze-out temperature Xf is typi- 
cally very small (xf = O(l/20)), justifying the iterative 
procedure Without further decay, the neutralino co- 
moving number density n x would have remained basi- 
cally constant. 

3. Neutralino decay into axino. In the scenario consid- 
ered here, the NLSP neutralino co-moving number den- 
sity after "freeze-out" will continue to decrease because 
of its decay P) into the axino LSP. This is presented in 
x < Xf and for T v <C H, n v is 
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Fig. § At Tj 
roughly given by 



n x (x) ~ n x eq (x f)C{x) exp 



dx' (r 



Xl % 



(4) 



where C{x) takes into account the temperature differ- 
ence between the photons and the decoupled neutralinos 
and (T x ) x is the thermally averaged decay rate for the 
neutralino at x, while H(m x ) = 2tt\/2ttNf/4:5 to! 2 /Mp. 
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FIG. 1. A schematic behavior of the co-moving number 
density: the thermal equilibrium (thick solid), NLSP neu- 
tralino (dash) and LSP axino (thin solid). 
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of the axino and the gaug- 
the neutralino is given by the 
term a Y C aY Y /{^V2^f a )[{^B a B a )ee + (** B & B & ) gg ] + 
otiCaWW I '(4v27T fa)[B — > W3], where $ is the chiral su- 
permultiplet containing the axion and the axino, while 
the vector multiplet B (W3) corresponds to the U(1)y 
(SU(2)l) gauge group with a coupling strength ay (02 )• 
The coefficients C a yy and C a ww are model dependent. 
Usually, one performs chiral transformations so that 
there is no axion- W IIV W ,1V coupling. This is equivalent 
to giving vanishing Peccei-Quinn charges to left-handed 
doublets. In this case, C a ww — and C a YY = C a 
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In 

the DFSZ model with (d c , e)-unification C a YY = 8/3, 
and in the KSVZ model for tQ = 0,-1/3, and 2/3 
CaYY — 0,2/3 and 8/3, respectively 17|. Below the 
QCD chiral symmetry breaking scale, C a77 and C a yy 
are reduced by 1.92. 

We first concentrate on the dominant decay channel (||) 
which is always allowed as long as m~ < m x . We will 



comment on other channels below, 
the process @ is given by 
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The decay rate for 
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where a em is the electromagnetic coupling strength, 
Caxi = {CaYY / cos 9w)Z\\, and N is a model depen- 
dent factor (N = 1(6) for the KSVZ (DFSZ) model). 

In the theoretically most favored case of a nearly pure 
B-ino fllglpl , the neutralino lifetime can be written as 



: 3.3 x 10" 2 sec 
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where the phase space factors from Eq. (|5j) have been 
neglected. 

A comment is in order regarding a plausible range 
°f fa PI- A somewhat model dependent lower bound 
f a > Kr GeV comes from astrophysical considerations, 
most notably from requiring that axions do not overly 
affect processes in globular clusters, red giants and in su- 
pernova 1987A. An upper limit of ~ 10 12 GeV is quoted 
in the context of cold axion energy density. A range 
10 9 " 10 GeV <f a < 10 12 GeV gives a cosmologically in- 
teresting values for the relic density of axions. 

Hard photons produced in the x decay thermalize 
via multiple scattering from background electrons and 
positrons (7 + e^7 + 7 + e) pp| , p^| . The process pro- 
ceeds rapidly for electromagnetic background tempera- 
tures above 1 MeV at which point background e — e pairs 
annihilate. To ensure efficient thermalization and in or- 
der to avoid problems with photodestruction of light ele- 
ments produced during nucleosynthesis, we require that 
the neutralino lifetime (H) is sufficiently less than about 1 
second (which also coincides with temperatures of about 
1 MeV). A modest requirement r < 10 _1 sec leads, in the 
case of the neutralino with a large B-ino component (a 
neutralino is never a pure B-ino state) , to an upper bound 
on f a which depends on m x . Also, at larger masses ad- 
ditional decay channels open up, most notably y — > Za. 
We can see that, for large enough values of the B-ino 
mass, the decay (||) will almost entirely take place before 
nucleosynthesis . 

The case of higgsino-dominated neutralino as the 
NLSP is probably less attractive and also more model 
dependent. First, the lifetime (j(^) is now be typically 
significantly larger, easily extending into the period of 
nucleosynthesis and beyond. This is caused by the sup- 
pression of the B-ino component through which the decay 
proceeds. (See the form of C axl below Eq. (jE))). Much 
lower values of f a could be considered as a remedy or 
much larger higgsino masses, and/or additional (model 
dependent) decay channels involving Higgs in the final 
state. 

In the MSSM, the higgsino relic abundance in the mass 
range allowed by LEP is typically very small, thus lead- 
ing to even smaller Sl~h 2 . One possibility would be to 
consider rather obese higgsino masses, above roughly 
500 GeV, where tt x h 2 > 1 again. A resulting value of 
il~h 2 would then depend on the actual size of the hig- 
gsino component in the decaying neutralino, as well as on 
the axion/axino model which would determine the cou- 
plings of the decay channels to the Higgs final state. One 
could also allow for a Higgs singlet and assume that its 
fermionic partner is mostly the NLSP. 

The resulting axino relic abundance today is simply 
given by 



a m x x 



(7) 
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since all the neutralinos have decayed into axinos. The 
axino will normally be produced relativistic, except when 
the ratio of the neutralino-axino mass difference to the 
axino mass is small, but will later redshift due to the ex- 
pansion of the Universe and become cold by the time of 
matter dominance. It is worth noting that the neutrali- 
nos will not dominate the energy density of the Universe 
before decaying; to see this we have to compare its life- 
time, given by Eq. (Q), with the time when the equality 
p x = p re i takes place. This time is easily computed ne- 
glecting the decay and amounts to 10 6 — 10 7 sec; we see 
therefore that the neutralinos never dominate the energy 
density and matter domination starts only after the pro- 
duced axinos become non-relativistic. 

With the lifetime (||) significantly larger than 10 _6 sec 
the neutralinos escape the high energy detectors. Thus 
phenomenology remains basically unchanged from the 
usual case where the neutralino is the LSP. In particu- 
lar, accelerator mass bounds on supersymmetric particles 
apply. But cases previously excluded by the constraint 
fl x h 2 < 1 can now be allowed via Eq. (@). This leads 
to a possibly dramatic relaxation of the parameter space 
of SUSY models. For example, in the MSSM the region 
of large higgsino masses mentioned above has been con- 
sidered cosmologically excluded but now can be again al- 
lowed if one takes a sufficiently small ratio m~/ m x . In the 
gaugino region it is normally reasonable to expect that, in 
order to reduce the LSP relic abundance below one |l8) , 
at least one sfermion with mass roughly below 500 GeV 
should exist. In the CMSSM, the same requirement of- 
ten provides an upper bound on the common scalar and 
gaugino masses of order one TeV over a large range 
of parameters. Both bounds which hold for a gaugino- 
like LSP away from annihilation resonances can now be 
readily relaxed. 

So far we have considered the neutralino as the NLSP. 
This assumption can easily be relaxed to accommodate 
any other superpartner, either neutral or carrying an 
electric or color charge, provided that its effective cou- 
pling with the axino is of order ~ l// a - All one needs 
to require is that the NLSP decay into the axino and 
the accompanying ordinary particle thcrmalization takes 
place before nucleosynthesis. If this can be achieved then 
cases previously considered excluded as corresponding to 
non-neutral LSPs can now again be allowed. 

In conclusion, we have shown that the axino can easily 
be the LSP with a mass in the GeV range. Such an 
axino would be a cold dark matter candidate for a natural 
range of the Peccei-Quinn scale f a . It is not impossible 
that, with or without its non-supersymmetric partner, 
the axino could dominate the matter in the Universe. 
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